Neutron star composition in strong magnetic fields 
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Abstract 

We study the neutron star composition in the presence of a strong magnetic field. 
The effects of the anomalous magnetic moments of both nucleons and electrons are 
investigated in relativistic mean field calculations for a /3-equilibrium system. Since 
neutrons are fully spin polarized in a large field, generally speaking, the proton 
fraction can never exceed the field free case. An extremely strong magnetic field 
may lead to a pure neutron matter instead of a proton-rich matter. 
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The discovery of pulsars [1] substantially stimulated the study of neutron stars, a 
class of extreme objects in astrophysics. Indeed, shortly after their discovery pulsars were 
identified as rotating magnetized neutron stars [2] with a surface magnetic field of 10 12 
- 10 14 G [3]. Recent observations of soft gamma repeaters (SGRs) have confirmed that 
they are newly born neutron stars with very large surface magnetic fields (up to 10 15 G) 
[4]. Such stars are named as magnetars [5]. It is expected that appreciably higher fields 
could exist in the interior of neutron stars although the total strength of internal magnetic 
fields remains unknown. The allowed internal field strength of a star is constrained by the 
scalar virial theorem [6, 7]. According to it, the maximum interior field strength could 
reach ~ 10 18 G. Naturally, one may expect considerable influences of such high fields on 
the properties of neutron stars. 

Theoretical investigation of a free electron gas in intense fields related to the neutron 
star crusts [7] and white dwarfs [8] as well as ideal noninteracting neutron-proton-electron 
(n-p-e) gas for the internal properties of neutron stars [9] has been carried out by sev- 
eral authors. The equation of state (EOS) of pure neutron matter under large magnetic 
fields was investigated in Ref. [10] where the strong interaction was taken into account. 
Theoretical calculations of spin polarized isospin asymmetric nuclear matter based on 
the Bruckner-Hartree-Fock formalism has been performed in Ref. [11]. Analogous cal- 
culations employing various mean-field models can be seen in Ref. [12]. The effects of 
strong magnetic fields on a dense n-p-e system under the beta equilibrium and the charge 
neutrality conditions were studied in Ref. [13]. Extremely large fields up to the proton 
critical field (here the critical field is defined by the condition that the particle's cyclotron 
energy is comparable to its rest mass [14]) B p c ~ 10 20 G and beyond were considered in 
their calculations. It was found that the nuclear matter in beta equilibrium practically 
converts into stable proton-rich matter. However, in these investigations the effects of the 
anomalous magnetic moments (AMM) were neglected, which may play an important role 
in the presence of such high fields. Recent studies demonstrated [15] that the contribu- 
tions from the AMM of nucleons become significant when the magnetic field B > 10 5 B^. 
Here B e c = 4.414 x 10 13 G is the electron critical field. Consequently, the softening of 
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the EOS caused by Landau quantization [13] is overwhelmed by the stiffening due to the 
incorporation of the nucleon anomalous magnetic moments. But a dramatic increase of 
the proton fraction with the increase of the magnetic field remains. In the study of Ref. 
[15] the AMM of electrons was not taken into account. 

The anomalous magnetic moment of an electron [16] may not be negligible due to 
its relatively large Bohr magneton. Because of the charge neutrality of the system, a 
significant effect on electrons will in turn influence protons. Considering that the proton 
fraction is crucial in determining the direct URCA process which leads to the cooling of 
neutron stars [17, 18], in this work we examine the problem with the inclusion of the 
AMM not only for nucleons but also for electrons. One may argue that the electron self- 
energy may not change substantially in magnetic fields when high-order terms are taken 
into account. However, systematic incorporation of high-order contributions beyond the 
AMM is not yet clear, which will be the topic of our forthcoming works. Here the effects 
of magnetic fields on different particles within the same system are taken into account at 
the equal footing. We consider the magnetic fields up to 1Q 6 B^. Though we do not yet 
know any object generating such fields, some possibilities have been suggested [19]. 

For the n-p-e system in a uniform magnetic field B along the z axis, the interacting 
Lagrangian can be written as [20] 

C = - e^p 7 ^ - ^c 6 //^ 

- e^A» - ^K ef i B a^F^ - m e ]4> e + ^ad^a - U(a) 
- l -u^ + \mlu^ - ^R M „ • RT + imjR M • R", (1) 

where the conventional notation has been adopted. U(a) is the self-interaction part of 
the scalar field [21] 

U(a) = \mlo* + l -b{ ga af + \c{ ga af. (2) 

= (0, 0, Bx, 0) refers to a constant external magnetic field and a^ u = |[7 M , 7"]. /j-n and 
fiB are the nuclear magneton of nucleons and Bohr magneton of electrons; k p = 3.5856, 
n n = —3.8263 and K e = a/n are the coefficients of the AMM for protons, neutrons and 
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electrons, respectively. The third set of parameters in Table II of Ref. [22] is used as 
the nucleon coupling strengths. It turns out g a = 8.7818, g w = 8.7116, g p = 8.4635, 
bgl = 27.9060, cg A a = -14.3989. This yields a binding energy B/A = -16.3 MeV, 
saturation density po — 0.153 fm -3 and bulk symmetry energy a sym = 32.5 MeV. 

The general solutions for the Dirac equation with the inclusion of the AMM are in 
the form of ip{X) oc e~ lEt+iPyV+ipzZ (f) p (py,p z , x) for protons and electrons and ip{X) oc 
e~ %Et+%p ' x <f> n (p) for neutrons. We have derived the concrete expressions of the Dirac spinors 
4>p{p y ,Pz, x) and n (p) in the chiral representation [23]. The chemical potentials of protons 
and neutrons are defined as 
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P P = e p f + g^uo + -gpRo,o, 

Pn = £f + guj^O — -j&pRsfl- 

They are related to the respective Fermi momenta via following equations: 

2 = (3) 2 " (Vm* 2 + 2eBu + SA)\ (5) 

fe) 2 = (^/) 2 -K + ^) 2 . (6) 

Here A = — ^K b p N B; S — ±1 for spin-up and spin-down particles, v is the quantum 
number of Landau levels for charged particles [14]. The effective nucleon mass m* = 
Mjy — go-c. Here a, and i? ,o are the mean values of the scalar field, the time-like 
component of the vector field and the time-like isospin 3-component of the vector-isovector 
field in neutron star matter, respectively. They are obtained by solving the non-linear 
equations of the meson fields. The total scalar density and baryon density are ps = Ps+Ps 
and p = po + Po respectively, where 
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The summation of v runs up to the largest integer for which [k^^ gj is positive. For spin- 
up protons v starts from 1 while for spin-down protons 0. It should be pointed out that 
here the so-called spin up and spin down are just relative notes since the wave functions 
are no longer eigenfunctions of 3-component spin operator [23], mainly attributed to the 
coupling of the spin to the magnetic field. The electrons are assumed to move freely in 
the strong magnetic field. The chemical potential of electrons p e = e^. Its relation to 
the electron Fermi momentum as well as the definition of the electron density pg are the 
same as given in Eqs. (5) and (9) except that the corresponding quantities are replaced 
by the electron ones. Numerical calculations are performed under the constraints of the 
charge neutrality pg = Po an d the /9-equilibrium p n = p p + p e [24]. These two constraint 
equations together with three meson equations are solved self-consistently in an iterative 
procedure. 

Figure 1 displays the effective nucleon mass and proton fraction as functions of the 
baryon density and magnetic field strength. Two situations with and without taking 
into account the anomalous magnetic moments are distinguished. In agreement with the 
finding of Refs. [13, 15], the effective mass remains unaltered from the field free case 
when B < lO 5 !^. If the field is increased to 10 6 B^, the effects of magnetic fields cause 
the m*/M]\r to increase when the AMM is included (n ^ 0) and decrease when the AMM 
is neglected (k; = 0). The proton fraction is shown to increase with the increasing of 
the magnetic field strength in the case of k — 0. A considerable enhancement of Y p is 
observed at B = 10 6 B^. For even stronger magnetic fields, one may speak about proton- 
rich matter [13]. The situation, however, becomes completely different if the anomalous 
magnetic moments are taken into account. At B = 10 5 B^, Y p is a little enhanced at lower 
density and then quickly approaches the field free case as the density increases. In the 
case of B = 10 6 B^, no protons survive at p < 4p . Evident suppression of the proton 
fraction is exhibited for the typical density range of neutron stars. 

In order to understand the above results, in the upper panel of Fig. 2 we depict the 
proton fraction as a function of the magnetic field strength, where the effects of the 
AMM are neglected. It can be found that for a fixed density the proton fraction begins 
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to increase rapidly at a certain value of the magnetic field. This is the critical point 
where both protons and electrons are completely spin polarized (i.e., only one Landau 
level is occupied.) due to Landau quantization while no direct effects from the magnetic 
field act on neutrons. The Fermi energies of electrons and protons fall drastically. As a 
consequence, a large amount of neutrons converts to protons. The system approaches to 
a well proton-rich matter. The variation of the critical field with the baryon density is 
displayed in Fig. 3 as the dashed line. 

In the lower panel of Fig. 2 we show the results including the effects of the AMM. 
Different curves are related to the different baryon densities. For each curve there exist 
two turning points. The first one corresponds to the critical field of electron polarization, 
its variation with density is depicted as the dotted line in Fig. 3. As soon as the electron is 
fully polarized, its Fermi energy will be considerably reduced. This causes the conversion 
of neutrons to protons in a beta-equilibrium system. The proton fraction is therefore 
enhanced substantially. Because of the different coefficients for the AMM of electrons 
and protons, the complete spin polarization of protons is reached later than electrons as 
shown in Fig. 3 with the dash-dotted line. There are no evident signals for the proton 
spin polarization in the lower panel of Fig. 2 since the rapid increase of the proton fraction 
has been induced by the electron polarization. When the field is further increased, the 
neutron can be fully spin polarized because of the coupling of the neutron spin to the 
magnetic field, an effect totally caused by the anomalous magnetic moment. The inverse 
process, p — > n, starts at the second turning points appeared on the curves plotted in 
the lower panel of Fig. 2. This causes the Y p to drop rapidly with the increase of B. A 
pure neutron matter becomes possible at very large field strength. The critical field for 
neutron polarization as a function of the baryon density is given in Fig. 3 as the solid 
line. One can see that the neutron spin polarization happens at B > 10 5 B^. We have also 
considered the situation where the electron AMM is switched off but the nucleon AMM 
is taken into account. No enhancement of the proton fraction is observed compared to 
the case that both the electron and nucleon AMM are neglected. While comparing with 
the condition that the AMM of both electrons and nucleons are considered, the proton 
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fraction is less suppressed. For instance, protons start to appear at p pa 2p rather than 
p pa 4p for B = 10 6 -B^. In order to check our numerical process we have revealed the 
results of Ref. [15], which is indicated in the lower panel of Fig. 1 as the long-dashed line. 
The above calculations have been repeated with other mean-field parameter sets used in 
Refs. [13, 15], the qualitative trend remains unchanged. 

In summary, we have studied the properties of neutron star matter under strong 
magnetic fields, with an emphasis on the neutron star composition. After taking into 
account the effects of the AMM of nucleons and electrons, the proton fraction is found to 
never exceed the field free case. Our results demonstrate that extremely strong magnetic 
fields may lead to a pure neutron matter rather than a proton-rich matter observed in 
Refs. [13, 15], mainly attributed to the complete neutron spin polarization induced by 
the AMM effects. 
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B=0, B=10 5 B^ 




Figure 1: The effective nucleoli mass m*/M N (upper panel) and proton fraction Y p (lower 
panel) as functions of the baryon density. Different curves are related to the different 
cases of magnetic fields B and with or without the inclusion of the anomalous magnetic 
moments k as indicated in the figure. In the lower panel, the long-dashed line reveals the 
case considered in Ref. [151. 

10 




Figure 2: The proton fraction as a function of the magnetic field strength. The upper 
panel shows the results without the AMM while in the lower panel the AMM is taken 
into account. 
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Figure 3: The critical field of particle spin polarization as a function of the baryon density. 
Different curves correspond to the different cases as indicated in the figure and discussed 
in the text. 
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